Empirical Bayes inference of pairwise F(ST) and its distribution in the genome

Populations often have very complex hierarchical structure. Therefore, it is crucial in genetic monitoring and conservation biology to have a reliable estimate of the pattern of population subdivision. F(ST)'s for pairs of sampled localities or subpopulations are crucial statistics for the exploratory analysis of population structures, such as cluster analysis and multidimensional scaling. However, the estimation of F(ST) is not precise enough to reliably estimate the population structure and the extent of heterogeneity. This article proposes an empirical Bayes procedure to estimate locus-specific pairwise F(ST)'s. The posterior mean of the pairwise F(ST) can be interpreted as a shrinkage estimator, which reduces the variance of conventional estimators largely at the expense of a small bias. The global F(ST) of a population generally varies among loci in the genome. Our maximum-likelihood estimates of global F(ST)'s can be used as sufficient statistics to estimate the distribution of F(ST) in the genome. We demonstrate the efficacy and robustness of our model by simulation and by an analysis of the microsatellite allele frequencies of the Pacific herring. The heterogeneity of the global F(ST) in the genome is discussed on the basis of the estimated distribution of the global F(ST) for the herring and examples of human single nucleotide polymorphisms (SNPs).     

Positive selection of Toll-like receptor 2 polymorphisms in two closely related old world monkey species, rhesus and Japanese macaques

Toll-like receptor 2 (TLR2) plays an important role in the recognition of a variety of pathogenic microbes. In the present study, we compared polymorphisms of TLR2 locus in two closely related old world monkey species, rhesus monkey (Macaca mulatta) and Japanese monkey (Macaca fuscata). By nucleotide sequencing of the third exon of TLR2 gene from 21 to 35 respective individuals, we could assign 17 haplotype combinations of 17 coding SNPs of ten non-synonymous and seven synonymous substitutions. A non-synonymous substitution at codon position 326 appeared to be differentially fixed in each species, asparagine for M. mulatta whereas tyrosine for M. fuscata, and may contribute to certain functional properties because it locates in the region contributing to ligand binding and interaction with dimerization partner of TLR2-TLR1 heterodimeric complex. Although TLR2 alleles have diverged to similar extent in both species, they have evolved in significantly different ways; TLR2 of M. fuscata has undergone purifying selection while the membrane-proximal part of the extracellular domain of M. mulatta TLR2 exhibits higher rates of non-synonymous substitutions, indicating a trace of Darwinian positive selection.     

Evolution of sex chromosomes ZW of Schistosoma mansoni inferred from chromosome paint and BAC mapping analyses

Chromosomes of schistosome parasites among digenetic flukes have a unique evolution because they exhibit the sex chromosomes ZW, which are not found in the other groups of flukes that are hermaphrodites. We conducted molecular cytogenetic analyses for investigating the sex chromosome evolution using chromosome paint analysis and BAC clones mapping. To carry this out, we developed a technique for making paint probes of genomic DNA from a single scraped chromosome segment using a chromosome microdissection system, and a FISH mapping technique for BAC clones. Paint probes clearly identified each of the 8 pairs of chromosomes by a different fluorochrome color. Combination analysis of chromosome paint analysis with Z/W probes and chromosome mapping with 93 BAC clones revealed that the W chromosome of Schistosoma mansoni has evolved by at least four inversion events and heterochromatinization. Nine of 93 BAC clones hybridized with both the Z and W chromosomes, but the locations were different between Z and W chromosomes. The homologous regions were estimated to have moved from the original Z chromosome to the differentiated W chromosome by three inversions events that occurred before W heterohcromatinization. An inversion that was observed in the heterochromatic region of the W chromosome likely occurred after W heterochromatinization. These inversions and heterochromatinization are hypothesized to be the key factors that promoted the evolution of the W chromosome of S. mansoni.     

Characterization of Streptococcus criceti insertion sequence ISScr1

A novel insertion sequence element of the IS982 family, ISScr1, was previously identified in Streptococcus criceti strain E49 as a disrupted paaB gene encoding an antigen I/II homologous protein. In this study, we identified two divergent inserted regions of ISScr1 in S. criceti E49 by inverse polymerase chain reaction, the gene-walking method, and screening of the partial plasmid library. Nucleotide sequence analysis indicated the possible explanation that transposition generated 8- and 9-bp direct repeat sequences. DNA hybridization analysis revealed that an identical hybridization pattern to ISScr1 was observed in the four S. criceti strains studied and that at least three copies of ISScr1 were preserved in S. criceti strains. In addition, we found different susceptibility to erythromycin and diverse agglutination properties induced by dextran in S. criceti. Furthermore, DNA hybridization analysis showed that no ISScr1-like copy was detected in the other 14 strains of oral streptococci tested.     

A novel composite retrotransposon derived from or generated independently of the SVA (SINE/VNTR/Alu) transposon has undergone proliferation in gibbon genomes

The superfamily Hominoidea (hominoids) comprises two families: Hominidae (hominids) and Hylobatidae (gibbons, also called small apes). The SVA transposon is a composite retrotransposon that occurs widely in hominoids and is considered to have been generated by stepwise fusions of three genetic elements: SINE-R, a variable number of tandem repeat (VNTR) sequence, and Alu. We identified a novel transposon whose basic structure is the same as that of SVA, with one prominent difference being the presence of part of prostaglandin reductase 2 (PTGR2) in place of SINE-R. We designate this composite transposon as PVA and propose two possible mechanisms regarding its generation. One is the derivation of PVA from SVA: the SINE-R region of SVA was replaced with a PTGR2 fragment by template switching. The other is the formation of PVA independently of SVA: a PTGR2 fragment was fused to an evolutionary intermediate comprising the VNTR and Alu regions. The nucleotide sequence of the junction between the VNTR and PTGR2 regions supports the second hypothesis. We identified PVA in the white-cheeked gibbon Nomascus leucogenys by analysis of genome sequence databases, and subsequent experimental analysis revealed its presence in all four gibbon genera. The white-cheeked gibbon harbors at least 93 PVA copies in its haploid genome. Another SVA-like composite transposon carrying parts of the LINE1 and Alu transposons in place of SINE-R, designated as LAVA, has recently been reported. The significance of the discovery of PVA is that its substituted fragment originates not from a transposon but from a single-copy gene. PVA should provide additional insights into the transposition mechanism of this type of composite transposon; the transposition activity is conferred even if the substituted fragment is not related to a transposon.     

Preferential paternal origin of microdeletions caused by prezygotic chromosome or chromatid rearrangements in Sotos syndrome

Sotos syndrome (SoS) is characterized by pre- and postnatal overgrowth with advanced bone age; a dysmorphic face with macrocephaly and pointed chin; large hands and feet; mental retardation; and possible susceptibility to tumors. It has been shown that the major cause of SoS is haploinsufficiency of the NSD1 gene at 5q35, because the majority of patients had either a common microdeletion including NSD1 or a truncated type of point mutation in NSD1. In the present study, we traced the parental origin of the microdeletions in 26 patients with SoS by the use of 16 microsatellite markers at or flanking the commonly deleted region. Deletions in 18 of the 20 informative cases occurred in the paternally derived chromosome 5, whereas those in the maternally derived chromosome were found in only two cases. Haplotyping analysis of the marker loci revealed that the paternal deletion in five of seven informative cases and the maternal deletion in one case arose through an intrachromosomal rearrangement, and two other cases of the paternal deletion involved an interchromosomal event, suggesting that the common microdeletion observed in SoS did not occur through a uniform mechanism but preferentially arose prezygotically.     

TNF-alpha and IL-4 regulate expression of IL-13 receptor alpha2 on human fibroblasts

Two interleukin 13 receptors (IL-13Rs) have been identified as IL-13Ralpha1 and IL-13Ralpha2. IL-13Ralpha1 is composed of a heterodimer consisting of IL-13Ralpha1 and IL-4 receptor alpha (IL-4Ralpha) as a signaling subunit. In contrast, IL-13Ralpha2 is known as a decoy receptor for IL-13. In this study, we investigated the expression of IL-13Rs on human fibroblasts. IL-13Ralpha2 was significantly up-regulated after stimulation with tumor necrosis factor-alpha (TNF-alpha) and/or IL-4. In contrast, IL-13Ralpha1 was constitutively detectable and was not up-regulated. After the induction of IL-13alpha2 by IL-4, STAT6 phosphorylation through IL-13Ralpha1 by IL-13 was inhibited. We also detected large intracellular pools of IL-13Ralpha2 in fibroblasts quantitatively. Furthermore, mobilization of the IL-13Ralpha2 protein stores from the cytoplasm to the cell surface was prevented by an inhibitor of protein transport, brefeldin-A. These results indicate that TNF-alpha and IL-4 synergistically up-regulate the expression of IL-13Ralpha2 decoy receptor on human fibroblasts by inducing gene expression and mobilizing intracellular receptors, and thus may down-regulate the IL-13 signaling.     

Impact of oxidative stress in aged mouse oocytes on calcium oscillations at fertilization

In vivo post-ovulatory aging of oocytes significantly affects the development of oocytes and embryos. Also, oocyte aging alters the regulation of the intracellular calcium concentration, thus affecting Ca(2+) oscillations in fertilized oocytes. Because reactive oxygen species (ROS) are known to significantly perturb Ca(2+) homeostasis mainly through direct effects on the machinery involved in intracellular Ca(2+) storage, we hypothesized that the poor development of aged oocytes that may have been exposed to oxidative stress for a prolonged time might arise from impaired Ca(2+)-oscillation-dependent signaling. The fertilization rates of aged oocytes and of fresh oocytes treated with 100 microM hydrogen peroxide (H(2)O(2)) for 10 min were significantly lower than that of fresh oocytes. Comparing within the fertilized oocytes, blastocyst formation was decreased while embryo fragmentation was increased similarly in the aged and H(2)O(2)-treated fresh oocytes. The frequency of Ca(2+) oscillations was significantly increased whereas the amplitude of individual Ca(2+) transients was lowered in the aged and H(2)O(2)-treated fresh oocytes. The rates of rise and decline in individual Ca(2+) transients were decreased in these oocytes, indicating impaired Ca(2+) handling. When lipid peroxidation was assessed using 4,4-difluoro-5-(4-phenyl-1,3-buttadienyl)-4-bora-3a, 4a-diaza-s-indacene-3-undecanoic acid (C11-BODIPY) in unfertilized oocytes placed in a 5% CO(2) in air atmosphere, the green fluorescence (indicating lipid peroxidation) increased faster in the aged oocytes than in the fresh oocytes. Furthermore, the green fluorescence in the aged oocytes was already approximately 20 times higher than that in the fresh oocytes at the beginning of the measurements. These findings support the idea that Ca(2+) oscillations play a key role in the development of fertilized aged oocytes.